Quantum transport in a normal metal/odd-frequency superconductor junction 
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Recent experimental results indicate the possible realization of a bulk odd-frequency superconducting state 
in the compounds CeCu2Si2, and CeRhlns. Motivated by this, we present a study of the quantum transport 
properties of a normal metal/odd-frequency superconductor junctions in a search for probes to unveil the odd- 
frequency symmetry. From the Eliashberg equations, we perform a quasiclassical approximation to account 
for the transport formalism of an odd-frequency superconductor with the Keldysh formalism. Specifically, we 
consider the tunneling charge conductance and tunneling thermal conductance. We find qualitatively distinct 
behaviour in the odd-frequency case as compared to the conventional even-frequency case, in both the electrical 
and thermal current. This serves as a useful tool to identify the possible existence of a bulk odd-frequency 
superconducting state. 

PACS numbers: 74.45.+C, 74.20.Rp,74.50.+r 



I. INTRODUCTION 

The symmetries of the superconducting order parameter 
with respect to orbital-, time-, and spin-space are governed 
by the Pauli principle. A wavefunction describing two elec- 
tronic states must be totally antisymmetric under exchange 
of the particle-coordinates. This leads to a finite number of 
allowed combinations for the symmetries of the wavefunc- 
tion. In a wide variety of superconductors ranging from those 
described with Bardeen-Cooper-Schrieffer/Eliashberg theory 
via spin-triplet superconductivity in ■^He, to strong-coupling 
superconductivity in high-Tc cuprates, the wave function of 
Cooper-pairs is even in the frequency domain. For such 
even-frequency pairing, the wavefunction may be even or odd 
in space depending on whether the Cooper-pairs form spin- 
singlets or -triplets. However, more exotic types of pairings 
than what is found in this wide range of materials, are in prin- 
ciple permitted. 

Recently, it was predicted that in a feiTomag- 
net/superconductor structure, a so-called odd-frequency 
pairing could take placai. Thus, the Cooper pair wavefunction 
is symmetric under exchange of spatial- and spin-coordinates, 
but antisymmetric under exchange of time-coordinates. This 
state had been proposed to exist by Berezinskii- a few decades 
earlier in the context of liquid '^He, and strong experimental 
evidence for odd-frequency pairing now exists^*. The study 
of such pairing in feiTomagnet/conventional superconductor 
junctions has been addressed by a number of authors over the 
last years'*. Furthermore, it was very recently predicted^'^'^ 
that due to spatial variation of the pair potential near a 
normal/superconductor (N/S) junction, the odd-frequency 
pairing state can be induced even in a conventional ballistic 
N/S system without spin-triplet ordering. The generation of 
different symmetry components and their effect on electrical 
transport in a normal/superconductor interface has also been 
studied in the diffusive limit'^ in the context of the proximity 
effect in unconventional superconductors'^ '". 

An issue that arises in the context of the odd-frequency 
pairing state, is if it can be realized in a bulk superconduc- 



tor, i.e. without a proximity effect. There have been several 
theoretical proposals for this in strongly correlated systems 
up to now"''^. To explore an odd-frequency pairing state in 
heavy-fermion superconductors is an interesting topic, and an 
assessment of the experimental properties of CeCu2Si2, and 
CeRhIn5 concluded that odd frequency pairing may be real- 
ized in these heavy-fermion compounds^. 
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FIG. 1: (Color online) Overview of the different symmetry states 
we will consider in the superconducting part of the clean, two- 
dimensional normal/superconductor junction. 

However, only a very limited amount of studies have ad- 
dressed the issue of identifying the odd-frequency pairing 
state in a bulk superconductor so fa r''''"*''^ . Hence, further 
clear-cut predictions are needed. 

In this paper, we present the quantum transport proper- 
ties of a normal metal/odd-frequency superconductor junction 
in the clean limit. We calculate the electrical and thermal 
conductances within the Blonder- Tinkham-KIapwijk (BTK) 
frameworbi^ taking account of the anisotropy of the pair 
potentialii. Our starting point is the Eliashberg equations that 
take into account the frequency-dependence of the pair po- 
tential. This constitutes a wide range of experimental predic- 
tions, which are routinely used to characterize superconduct- 



ing statei 
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Our main result is that the odd-frequency 



symmetry affects the charge (thermal) transport in an essen- 
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tial manner at low energies (temperatures). This provides a 
useful tool in identifying this highly unusual superconducting 
state. 

To elucidate the physics in a transparent manner, we em- 
ploy a simple two-dimensional calculation in the clean limit. 
We approximate the superconducting gap with a step-function 
in space, which in the isotropic even-parity s-wave case 
should be an excellent approximation for low-transmission 
barriers. Since the low-transmission case probably is the 
most realistic scenario experimentally, we restrict our atten- 
tion to this. In the anisotropic even-parity and odd-parity 
cases (corresponding e.g. to the high-T^ superconductors and 
Sr2Ru04), the gap may undergo a severe depletion near the 
barrier even for low-transmission interfaces due to the forma- 
tion of zero-energy states^^. The method used in this paper 
may still be able to capture qualitative features of the trans- 
port properties even in those cases, just as in the case of the 
d-wave superconductors'^. Our results are in fact consistent 
with recent findings'' including a self-consistent solution of 
the spatial variation of the superconducting gap near the inter- 
face. 

We will use boldface notation for 3-vectors, . ? . for 4x4 
matrices, and for 2x2 matrices. Pauli-matrices in particle- 
hole x spin (Nambu) space are denoted as pi, while Pauli- 
matrices in spin-space are written as r^. 



11. THEORETICAL FORMULATION 

A. Equations for odd-frequency superconductivity 

The frequency-dependence of the superconducting order 
parameter may be naturally taken into account in the approach 
developed by Eliashberg^^, where details of the electron- 
boson interaction are taken seriously. This contrasts the usual 
weak-coupling picture where the pairing interaction is taken 
to be constant. For our purposes, the following Hamiltonian 
is an appropriate starting point: 

^^E/ dr^a(r)^i?/(r)Va(r)+ j drbHr)Hb{r)b{r) 

a. 

+ Y.j /"drdr'V(r-r')Vl(r)V'a(r)[6(r)+6t(r)], (1) 

a 

where Hf is the Hamiltonian for free fermions which we as- 
sume may be written as Hf{r) ~ ^2m(^ ^ leA)^ — p, 
while Hb is the Hamiltonian for free bosons. Above, a de- 
notes the spin index while and b are fermion and boson op- 
erators, respectively. Introducing the Fourier-transformation 
Kr) = ^Eq^e-''^■^ = b^ + bl^, we obtain the 
Heisenberg equations of motion 

idtipa{r,t) = iJ/(r)-0„(r,t) +^C(i",t,q)V'a(r), 

q 

idtijiir,t) = -i/;(r)^t(r,t) _^^(r,i,q)V,t(r), (2) 

q 



where (^{r,t,q) = Vqi3q(t)e~"i '" and Vq is the Fourier- 
transform of V(r — r'). Note that V is not the effective pair- 
ing potential between electrons. Having obtained the time- 
derivatives of the fermion operators, we may now calculate the 
equation of motion for the Green's functions. This procedure 
is standard and covered in e.e. Refs. l24ll25ll26ll27i Taking into 
account the effect of the electron-boson interactions explicitly 
in the Hamiltonian naturally includes a frequency-dependence 
in the effective electron-electron interaction^^ which is ob- 
tained by integrating out the bosonic degrees of freedom in 
the partition function. The effective electron-electron interac- 
tion mediated by a boson excitation may in general be written 
as 

n<,.a,.^. ,3, 

where q = k — k' and fl ~ oj — oj' are the momen- 
tum and energy transfers, respectively, of the interaction pro- 
cess. Above, cjq is the frequency of the boson propagator 
Note that the pairing potential in Eq. (O is even in il, i.e. 
V{q, il) — y(q, — fi). The self-consistency equation for the 
order parameter quite generally has the structuroii 

which may be re-written as 

by exploiting V{q,ft) = T^(q, — O). The above equa- 
tions show that both A(k, a;) = A(k, —uj) and A(k, — 
— A(k, —uj) are possible solutions of the gap equation. There- 
fore, although the pairing interaction is even in frequency, the 
gap A in principle may be both even or odd in frequency. 
In fact, it is in general a superposition of even- and odd- 
frequency components"'^**. Assuming that the energy trans- 
fer is small compared to the term containing the momenta in 
Eq. (O, \ ft\ ^ jtiiql, one obtains a part of the pairing poten- 
tial which is linear in ui and uj' and one that is quadratic in 
the same quantities" . The former part is the necessary in- 
gredient to obtain a superconducting order parameter that is 
odd-in-frequency. It is also possible to adopt a purely phe- 
nomenological approach to an odd-frequency superconductor 
by assuming the frequency-dependence of the gap a priori^. 

Let us now consider the structure of the Green's function 
matrix for an odd-frequency superconductor. It is instructive 
to briefly mention the result for an ordinary BCS superconduc- 
tor, which has an even frequency-singlet-even parity (ESE) 
symmetry. In the BCS case, one obtains 

(i J^P3 - e - A(ri)) G^{1, 2) - 6(1 - 2)1. (6) 

Assuming a homogeneous and isotropic system where the 
Green's function only depends on the relative coordinates 
t = ti — t2 and r = ri — r2, and where A(ri) = A is 
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a constant, one may Fourier-transform Eq. (|6]l according to 
GR(p, e) = / / dre-'P''dte'^*G'^(r, t), where e and p is the 
quasiparticle energy measured from Fermi level and momen- 
tum, respectively. We then obtain 

(£P3-ep-A)G'^(p,e) = i, (7) 

which upon matrix inversion yields the well-known BCS so- 
lution. The quasiclassical Green's functions g{pp,r;e,t) is 
obtained from the Gor'kov Green's functions G(p, r; e, t) by 
integrating out the dependence on kinetic energy, assuming 
that G is strongly peaked at Fermi level, 

g{pp, r;e,t)^ ^ J d^pG(p, r; e, t). (8) 

The above assumption is typically applicable to superconduct- 
ing systems where the characteristic length scale of the per- 
turbations present, namely superconducting coherence length, 
is much larger than the Fermi wavelength. The correspond- 
ing characteristic energies of such phenomena must be much 
smaller than the Fermi energy ep. The quasiclassical Green's 
functions may be divided into an advanced (A), retarded (R), 
and Keldysh (K) component, each of which has a 4 x 4 matrix 
structure in the combined particle-hole and spin space. One 
has that 

where the elements of .g(pF, r; e, t) read 

The quantities g and / are 2 x 2 spin matrices, with the struc- 
ture 

g^i^^^ . (11) 

Due to internal symmetry relations between these Green's 
functions, all of these quantities are not independent. In par- 
ticular, the tilde-operation is defined as 

/(pF,r;e,i) = /(-pF,r; -e,t)*. (12) 

For a bulk s-wave superconductor, the retarded part may be 
expressed in terms of the normal and anomalous Green's func- 
tions g and / as follows: 

^ Uil2C-'x -gl j' 

Here, x is the globally broken U(l) phase associated with 
the spontaneous symmetry breaking of the superconducting 
state. In the odd-frequency case, however, one finally arrives 
at 

[e/53-ep-A(£)]G'^(p,e) = i, (14) 



where now A(e) is the odd-frequency gap matrix. Note that 
Eq. ( fT4l l is equivalent to the well-known Eliashberg equation. 
The structure of the Green's function for an odd-frequency su- 
perconductor may be different from Eq. (fT3l i depending on the 
spin-symmetry. For instance, the bulk Green's function matrix 
for an odd-frequency spin-triplet even-parity superconductor 
has the structure: 

Performing a quasiclassical approximation on Eq. (fT4l i yields 
the Eilenberger equation, which reduces to the Usadel equa- 
tion in the dirty limit. Note that for both even- and odd- 
frequency superconducting order parameters, the pairing in- 
teraction itself is always even in the frequency coordinate. 

A quite general formalism for treating quantum trans- 
port in non-uniform superconducting systems, e.g. nor- 
mal/superconductor heterostructures, has been developed by 
Tanaka and co-workers'"'. For instance, the conductance 
spectra of a normal/superconduct or J unction may be ob- 
tained along the lines of Refs. [30ll3 ll by numerically solv- 
ing the Usadel equation using Nazarov's generalized bound- 
ary conditionsiS. Interestingly, taking the Umit Rd and 
6* ^ in this formalism, where Rd represents the resistance 
of the normal metal region and is a measure of the prox- 
imity effect, leads to the well-known expression for the con- 
ductance obtained in the BTK-formalism'^. This may be seen 
specifically for the electrical conductance by consulting Eqs. 
(15) and (16) in Ref. l30l and for the thermal conductance in 
Eq. (19) of Ref. [U Therefore, since the above treatment of 
the Eliashberg equation shows that the odd-frequency depen- 
dence of the gap may be taken into account simply by sub- 
stituting A A(e), quantum transport for an odd-frequency 
superconductor can be treated in the BTK-formalism by per- 
forming the same substitution. However, the derivation of 
the Bogolioubov-de Gennes equation for odd-frequency su- 
perconductivity is challenging since it is not obvious how to 
take into account the strong retardation effects of the pairing 
potential. 



B. Transport formalism 

We adopt the Keldysh formalism using Nazarov's general- 
ized boundary conditions '° to obtain the electrical and ther- 
mal conductance for odd-frequency superconductors. We as- 
sume, without loss of generality, that the gap A(e, i?) has an 
opposite-spin pairing symmetry in both the singlet and triplet 
case. To encompass accessible experimental techniques, we 
will focus on two experimentally accessible quantities that 
encode how the odd-frequency pairing symmetry is mani- 
fested in transport properties: namely, the normalized charge- 
conductance G{eV) for T = and the thermal-conductance 
k{T). The procedure for obtaining these quantities is treated 
in detail in Refs. I30ll3lll32t In the limit of zero resistance in 
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the normal part and vanishing proximity effect, one finds 



G 



1 

Gn 

7r/2 
-7r/2 



r/2 



7r/2 



d7?cost9r+(ey,z9), 



-2/32r_(eV,i?) 



dz9de- 



4Ao cosh"' (cos z?)- 



(16) 



where Gn is the normal-state conductance and we have de- 
fined 



r„(£,i?) = 1 + 



-17+ 



r2_e'(T--T+) 



(17) 



Above, we have introduced ?9+ = t?, ?9_ = tt — i?, and 
n± = ^(l^+sign(e)/g±)/2, where sign(e) -sign(e) 
for 0^0. The phase of the superconducting gap is 
contained in the factor e''''± = c''^(''±) = /±/|/±|. The 
quantities g± and f± are the asymptotic values of the nor- 
mal and anomalous Green's functions of the odd-frequency 
superconductor in a gauge where the superconducting or- 
der parameter is real: g± ~ e/y/s^ — |A(£,t?±)p, f± = 
A(e,'(?±)/v^|A(£, ??±)|2 - We have introduced = 
~\Z/ cos 1?, where Z denotes the strength of the scattering po- 
tential near the barrier. In what follows, we fix Z = 3, corre- 
sponding to a typical low-transparency barrier which is exper- 
imentally realistic. Note that in the expression for k, we have 
considered the linear response regime for a small temperature- 
gradient in the system and introduced j3 = 1/T where T is the 
temperature of the reservoirs. 
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FIG. 2: (Color online) Plot of G for isotropic, even parity super- 
conductors and odd parity superconductors for both even and odd 
frequency pairing. 



m. RESULTS AND DISCUSSION 

Depending on the symmetries with respect to sign inversion 
of frequency and momentum, corresponding to e ^ 
d ^ Ti + -d, the gap may be classified as seen in Tab. I] In 
each case, we will model the gap A(e, d, T) as illustrated in 
the table. In the angular dependence of the odd-parity gaps, 
a denotes the misorientation angle between the antinodes and 
the interface normal (see Fig. [T]). The motivation for mod- 
elling the frequency-dependence of the superconducting gap 



as we have done in Tab. Uis that it features the low-energy be- 
haviour of the proximity-induced odd-frequency gap in dirty 
ferromagnet/superconductor structures^-' and that it exhibits a 
similar energy dependence to the gap seen in strongly corre- 
lated electron systems considered in Ref.fisl 

Recently, it was demonstrated that the odd-frequency pair- 
ing is quite generally induced near the normal/superconductor 
interface by a fully self-consistent calculation of the supercon- 
ducting correlations^. In an ETO superconductor with a = 0, 
corresponding to perfect formation of zero-energy states, an 



5 



TABLE I: Overview of the specific gap forms we will consider in this 
paper. We model the temperature-dependence of A{T) with A{T) — 
Ao tanh(1.74^rc/r - 1), and = Ao/1.76. 



Symmetry Specific gap form A (e, i9, T) 
ESE A{T)A„ 
OTE A{T)e /[l + {e/Aof] 

OSO ^(r)ecos(i9-Q)/[l + (e/Ao)^] 
ETO A{T) Aq cosji^ -a) /[I + {e/Aof] 



OTE pairing is induced near the surface. Thus, the forma- 
tion of zero-energy states may be re-interpreted as a manifes- 
tation of the odd-frequency superconductivity near the inter- 



face. The odd-frequency symmetry may permit the existence 
of gapless single-particle excitations at Fermi level. On the 
other hand, when the nodal direction is parallel to the inter- 
face normal (a = tt/2), only the even-frequency states exist 
at the interface. 

In a similar manner, the OSO pairing state can be induced 
near the interface of a clean normal/superconductor junc- 
tion when the superconductor has an ESE symmetry. One 
may also apply this discussion to bulk odd-frequency super- 
conducting states. In this scenario, the ETO (ESE) pair- 
ing can be induced at the interface for an OTE (OSO) bulk 
superconductor^. This should have clearly observable conse- 
quences for the quantum transport properties of a normal/odd- 
frequency superconductor junction. We now proceed to inves- 
tigate this in further detail. 




FIG. 3: (Color online) Plot of k for isotropic, even parity super- 
conductors and odd parity superconductors for both even and odd 
frequency pairing. A power-law dependence with exponent ~ 1 is 
observed for both of the odd frequency symmetries (see insets). 



Consider first the left column of Fig. |2] where we have plot- 
ted G as a function of bias voltage for the even-parity sym- 
metries. In the even-frequency case, the usual singularity at 
eV ~ AqAq is present. In the odd-frequency case, we see a 
qualitatively different behaviour of the conductance. First of 
all, G at low bias voltage is greatly enhanced compared to the 
even-frequency case for < 1, and the formation of a zero- 
bias conductance peak is clearly seen at Aq ~ 1. For Aq > 1, 



the conductance is similar to the even-frequency case for a re- 
duced value of the gap. This may be understood as follows. 
For ^0 < 1, the inequality e > A(e) is satisfied for all e with 
our choice of gaps (Tab. |l|i. This corresponds to gapless su- 
perconductivity. For Aq > 1, the gap becomes larger than e 
below a certain (finite) value of e, similar to the BCS-gap. 

The middle and right columns of Fig. |2]shows G as a func- 
tion of bias voltage for the odd-parity symmetries. In the OSO 
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case, a gap like structure is seen at a = 0. This is because ESE 
pairing is induced near the interface due to the sign change of 
the pair potential.^ This ESE pairing is responsible for the gap 
like structure of the conductance spectra, similar to the ESE 
case in Fig. 2. In contrast, OSO pairing remains near the in- 
terface at a = 7r/2. Thus, a zero-bias conductance peak is 
seen. On the other hand, in the ETO case at a = 0, a zero- 
bias conductance peak appears due to the induced OTE pair- 
ing near the interface^, similar to the OTE case in Fig. 2. At 
a = 7r/2, ETO pairing survives near the interface and hence 
the even frequency character of the pair amplitude results in a 
V-like shape of the spectra. Interestingly, OSO and ETO cases 
have the opposite tendency although their t?-dependencies are 
the same. Furthermore, the sign change of the gap produces 
a qualitative difference in the spectra between OTE and OSO 
with a — Q junctions. Thus, G is phase sensitive not only 
in even frequency superconductor junctionsi^, but also in odd 
frequency superconductor j unctions. 

We next investigate the thermal conductance k, shown in 
Fig. [3] The left column corresponds to the even parity case, 
where the usual exponential dependence on T is recovered 
for the ESE case-^'*. In the OTE case, k mimics the ESE case 
for Aq > 1 just as for the charge conductance. Otherwise, 
power-law dependence with exponent ~ 1 is observed due to 
the node of the gap at zero energy. Thus, the nodes in the fre- 
quency domain of an isotropic odd-frequency superconduc- 
tor causes k to behave like in an anisotropic even-frequency 
superconductor. In the middle and right columns of Fig. |3] 
we give k in the odd parity case. The well-known result of 
exponential dependence for a = is recovered in the ETO 
case. The OSO case again displays power-law behaviour sim- 
ilar to the OTE case for Aq < 1. However, the exponential 
dependence again occurs for > 1 in the OSO case with 
a = 0. When a — tt/2, there is exclusively power-law de- 
pendence, with exponent ~ 1. While the OTE case only has 
nodes in energy, the OSO case has both nodes in energy and 
in k-space, but this does not appear to influence the exponent 
of the power-law dependence. 



IV. SUMMARY 

In summary, we have studied quantum transport in a nor- 
mal metal/superconductor junction, considering how a bulk 
odd frequency symmetry in the superconductor is manifested 
in the electrical and thermal conductance of the junction. The 
odd frequency symmetry is found to display qualitatively dis- 
tinct behaviour from the even frequency case. This reflects 
the fact that the electrical conductance is sensitive to the pres- 
ence of odd frequency pairing at the interface, whereas the low 
temperature behavior of the thermal conductance reflects the 
node of the gap in the frequency domain. Moreover, one may 
distinguish the even and odd parity cases for an odd frequency 
symmetry (OTE and OSO, respectively) by means of their dif- 
ferent characteristic tunneling spectra. Our predictions should 
be useful for a wide range of experimental techniques, and are 
thus a helpful tool in identifying the possible existence of bulk 
odd frequency superconductors, with CeCu2Si2, and CeRhlns 
currently presenting themselves as the most promising candi- 
dates. 
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